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EFFECT OF INTERACTION ON LANDING-GEAR BEHAVIOR AND DYNAMIC LOADS IN A 

FLEXIBLE AIRPLANE STRUCTURE 1 

By Fbancis E. Cook and Benjamin Mtlwitzky 


SUMMARY 

The effects of interaction between a landing gear and a flexible 
airplane structure on the behavior of the landing gear and the 
loads in the structure have been studied by treating the equations 
of motion of the airplane and the landing gear as a coupled 
system. The landing gear is considered to have nonlinear 
characteristics typical of conventional gears, namely, velocity- 
squared damping, polytropic air-compression springing, and 
exponential tire force-deflection characteristics. For the case 
where only two modes of the structure are considered, an equiva- 
lent three-mass system is derived for representing the airplane 
and landing-gear combination, which may be used to simulate 
the effects of structural flexibility in jig drop tests of landing 
gears. 

As examples to illustrate the effects of interaction, numerical 
calculations, based on the structural properties of two large air- 
planes having considerably different mass and flexibility charac- 
teristics, are presented. For the particular cases considered, 
it was found that the effects of interaction can result in appreci- 
able reductions in the magnitude of the landing-gear force, 
particularly when the flexibility of the airplane structure is large 
and the natural frequency is small. Thus, neglect of interaction 
effects, that is, the use of the landing-gear forcing function for a 
rigid airplane, in a dynamic analysis of a flexible airplane can 
lead' to the calculation of excessive loads in the airplane struc- 
ture. 

In the case of one of the airplanes considered, the structural 
loads calculated from the interaction solutions are greater 
than those for a completely rigid airplane treatment {rigid struc- 
ture subjected to rigid-body forcing function) because the effects 
of dynamic magnification more than overcome the reduction in 
landing-gear force due to interaction. In the case of the second 
airplane, because of the relatively large natural period of the 
structure in comparison with the duration of the impact pulse, 
the dynamic magnification factor is appreciably less than unity. 
This effect, coupled with the reductions in landing-gear force 
due to interaction, results in structural loads that are less than 
those for a rigid airplane. 

INTRODUCTION 

In the design of landing gears it is usually assumed that 
the airplane is a rigid body and development tests are 


frequently carried out in a drop-test jig with a landing gear 
attached to a concentrated mass. In so doing, it is tacitly 
assumed that the interaction between the motions of the 
landing gear and the deformations of the airplane structure 
has little or no effect on the behavior of the landing gear. 
Also, load time histories obtained on a rigid-body basis are 
often used as the forcing function in a dynamic analysis to 
determine the inertia loads and stresses in flexible airplane 
structures, again under the assumption that the behavior 
of the landing gear is independent of the effects of airplane 
flexibility. Although it has been recognized that this assump- 
tion is not altogether valid, the errors involved have not been 
considered particularly significant in the past because: 
(a) The errors were thought to be on the conservative side 
and (b) until comparatively recently main landing gears 
have generally been located very dose to the nodes of the 
fundamental bending mode of the wing, and the airplane 
therefore closely approximated a rigid body insofar as the 
behavior of the landing gear is concerned. However, the 
trend toward increased size of airplanes, the disposition of 
large concentrated masses in outboard locations in the 
wings, the use of thinner wings, and the development of 
unconventional configurations tend to increase the flexibility 
of the airplane structure and reduce the natural frequencies 
of vibration. These characteristics tend to cause an increase 
in the amplitudes of the oscillatory motions of the landing- 
gear attachment points relative to the center of gravity of 
the flexible system during impact so that the effects of inter- 
action are increased, both with regard to the behavior of the 
landing gear and the dynamic loads in the structure, par- 
ticularly when the natural period of the fundamental mode 
of the structure approaches the time duration of the impact 
pulse. 

A number of analytical studies and some simplified model 
tests (refs. 1 to 5) which have been made to evaluate the 
effects of structural flexibility on landing-gear loads have 
indicated some reduction in landing-gear force due to the 
effects of structural deformation. However, in view of the 
fact that these previous investigations considered only 
rather highly idealized linear-spring landing gears with 
either no damping at all or viscous damping, a further study 
of the effects of interaction between the landing gear and 


i Supersedes NAOA Technical Note 3467 by Francis E. Cook and Benjamin MUwitsky, 1955. 
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the airplane structure has been made with a more realistic 
representation of the landing gear. In the present analysis, 
as in reference 6, the landing gear is considered to have 
velocity-squared damping, polytjopic air-compression spring- 
ing, and exponential tire force-deflection characteristics, as 
is the case with conventional oleo-pneumatic landing gears 
in current use. The particular purposes of this investigation 
are to evaluate the effects of interaction on the behavior of 
the landing gear and to study the errors introduced into the 
calculated loads in the structure (applied loads, accelerations, 
bending moments, and shears) when a dynamic analysis is 
made on the basis of applying the landing-gear forcing func- 
tion for a rigid body to a flexible airplane. For these pur- 
poses, case-history studies, based on the structural properties 
of two large airplanes having considerably different mass and 
flexibility characteristics, are presented. In order to cover 
a range of parameters, the landing gear of each airplane was 
assumed to be located at three arbitrary spanwise positions 
in addition to its original location. Only symmetrical 
impact conditions are considered. 

The basic analysis of the landing gear and the airplane 
structure as a coupled system is presented in a general form. 
In the numerical examples presented, however, the system 
is simplified by considering the motions of the airplane in its 
first two structural modes only. With these restrictions, 
the combination of airplane and landing gear can also be 
represented by an equivalent three-mass system which may 
be used in jig drop tests of landing gears to simulate the 
primary effects of structural flexibility. A similar type of 
concentrated-mass system was used in the study of the 
hydrodynamic impact of a flexible seaplane in reference 7. 

SYMBOLS 

GENERAL 

g gravitational constant 

E l lift factor, L, 0 JW, 0t 

L, ol total lift force (half the airplane) 

t time after initial contact 

t time variable of integration 

T time to maximum landing-gear force 

t' time after maximum landing-gear force 

V Vg vertical velocity at initial contact 
W t0 1 total weight (half the airplane) 

circular frequency of sine pulse 
Hi circular frequency , of cosine pulse 

X any variable 

X p value of any variable X at end of pth interval sub- 

sequent to beginning of numerical integration 
procedure 

LANDING GEAR 

A a pneumatic area of shock strut 

A,, hydraulic area of shock strut, Ai-Aj, 

Ai internal cross-sectional area of shock-strut inner 

cylinder 

A n net orifice area of shock strut, Ao-A, 

A 0 area of fixed opening in orifice plate 
A p cross-sectional area of metering pin or rod in 

plane of orifice 


C a orifice discharge coefficient 

F Vf vertical component of force in shock strut subse- 

quent to beginning of shock-strut deflection 
Fr t (zu) vertical force applied to tire at ground 
m, unsprung mass below shock strut 

m,r,m' constants in tire force-deflection relationship 
n polytropic exponent for air-compression process in 

shock strut 

p„ o air pressure in shock strut when fully extended 
p mass density of hydraulic fluid 

v 0 air volume of shock strut when fully extended 

s shock-strut stroke 

t { duration of impact pulse 

0 angle between shock-strut axis and vertical 

W u weight of unsprung mass below shock strut 
Zf vertical displacement of landing-gear attachment 

point from position at initial contact 
z K vertical displacement of axle from position at initial 

contact 

DISTRIBUTED STRUCTURE 

a n generalized coordinate for nth mode 

<p angle of twist of transverse station 

ip„ modal function for torsion in nth mode 

w vertical displacement of elastic axis from position at 

initial contact 

w n modal function of elastic axis for bending in nth 
mode 

f vertical displacement of station mass centers from 

position at initial contact 

modal function of station mass centers for coupled 
bending torsion in nth mode 

£„ modal amplitude of landing-gear attachment point 

for coupled bending-torsion in nth mode 
e chordwise distance between elastic axis and station 

mass center 
b wing span 

M b bending moment 

F vertical component of applied landing-gear forco 

fi natural frequency of first deflection mode 

I 0 polar moment of inertia of wing cross section about 

station mass center 

/„ polar moment of inertia of wing cross section about 
elastic axis 

K radius of gyration of 'wing station about elastic axis 
L lift force per unit length of span 

m mass per unit length of span 

M n generalized mass for nth mode (half the airplane) 
w n circular frequency of nth mode 

Q n generalized force in nth mode 

S shear 

t n natural period of nth mode 

x chordwise distance between elastic axis and any 

arbitrary point 

x f chordwise distance between elastic axis and landing- 

gear attachment point 

y spanwise distance from airplane center plane to any 

transverse station 
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y g spanwise distance from, airplane center plane to 

landing-gear station 

z vertical displacement of any point from position at 

initial contact 

Zj vertical displacement of landing-gear attachment 

point from position at initial contact 
z u vertical displacement of axle from position at initial 

contact 

5a* virtual displacement of generalized coordinate of 
nth mode 

Sll 7 * virtual work in nth mode 

EQUIVALENT THREE-MASS SYSTEM 

a Q vertical displacement of center of gravity of spring- 

connected masses from position of initial contact 
k spring constant 

Lf lift force acting on mass m f 

L, lift force acting on mass m, 

m f mass acting directly on landing gear 
m, elastically supported mass 

wj natural frequency of vibration of spring-connected 

masses 

u deflection of spring 

W/ weight of mass acting directly on landing gear 

W, weight of elastically supported mass 

z f vertical deflection of landing-gear attachment point 

2 S vertical deflection of elastically supported mass 

z u vertical displacement of axle from position at initial 

contact 

AERODYNAMIC 

A total iving area (half the airplane) 

At wing area assumed concentrated at station i 
Cz, lift coefficient 

C ho lift coefficient at instant of initial contact 

Cjj a lift-curve slope 

7 flight-path angle 

7 „ flight-path angle at instant of initial contact 

p mass density of air 

Vl landing speed of airplane 

Subscripts: 
a aerodynamic 

/ landing-gear attachment point 

g landing-gear station 

i any spanwise station 

n pertaining to the nth mode 

0 zero or rigid-body mode 

r at instant of initial shock-strut motion 

T at instant of maximum landing-gear force 

max maximum 

The use of dots over symbols indicates differentiation with 
respect to time t or r. All translations are positive downward 
(see figs. 1 to 3). The absolute value of any term ( ) is 
indicated by |( )|. 

ANALYSIS 

In order to study the behavior of a landing gear and a 
flexible airplane structure os mutually interacting elements of 
a coupled system, the equations for the landing-gear 


force are combined with the equations of motion of the 
structure. The motions of the structure are treated by the 
mode-superposition approach, wherein the deflections of the 
structure are expanded in terms of its natural modes of 
vibration. The effects of interaction between the landing 
gear and the structure are introduced by expressing the 
landing-gear force in terms of the motions of the landing-gear 
attachment point and the wheel axle (or unsprung mass) 
rather than as an arbitrary function of time. 

Because conventional oleo-pneumatic shock struts do not 
begin to deflect until some finite time after initial contact of 
the tire with the ground, the impact is treated in two parts, 
namely, the phases prior to and subsequent to the beginning 
of shock-strut deflection, where the initial conditions for the 
second phase are determined from the terminal conditions 
for the first phase. 

In the first part of the analysis, the equations for the 
landing-gear force are presented. Then, the deflections of 
the structure are expanded in terms of coupled modes and the 
resulting equations of motion for the system are presented in 
a general form. For the purpose of indicating the quanti- 
tative effects of interaction, however, the system used in the 
numerical trend studies has been simplified by restricting 
consideration of the structural deflections to the first two 
modes of the expansion. Within the framework of this 
two-mode treatment, it is also shown that the airplane struc- 
ture can be represented by an equivalent system of spring- 
connected concentrated masses, which may be used to 
simulate the effects of structural flexibility in jig drop tests 
of landing gears. 

LANDING-GEAR FORCE 

An analysis of the behavior of the conventional type of 
oleo-pneumatic landing gear was presented in reference 6. 
In this study the mass above the landing gear was considered 
as a rigid body; the system treated therefore had two degrees 
of freedom and is schematically represented in figure 1 . The 
analysis of the landing gear considered the velocity-squared 
damping of the metering orifice, the air-compression springing 
of the shock strut, the nonlinear force-deflection charac- 
teristics of the tire, and the internal shock-strut friction 
forces. Calculated time, histories of the landing-gear forces 
and the motions of the system were in good agreement with 
experimental data obtained in drop tests. 

In the present study the rigid mass is replaced by a flexible 
airplane structure, but the treatment of the landing gear is 
essentially the same as that in reference 6. However, since 
conventional landing gears are inclined forward so as to 
minimize normal forces and bending moments due to the 
combination of vertical and drag forces, it will be assumed 
that the resultant force on the landing gear lies along the axis 
of the shock strut so that strut bending moments and 
resulting internal friction forces are neglected in the present 
analysis. 

In view of the fact that conventional oleo-pneumatic shock 
struts are preloaded with air and therefore do not begin to 
deflect until some finite time U after initial contact of the 
tire with the ground, the impact must be treated in two 
phases. In the first phase, since the strut is effectively rigid, 
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(a) System with two degrees of freedom. 

Fiquhe 1. — Dynamical system (rij 

the landing gear has only one degree of freedom and the 
motion of the complete system of the landing gear and 
airplane is governed by the force between the tire and the 
ground. This ground force arises from the deflection of the 
tire and, in general, may be written as 

Fv t =Fv t (z„) ( 1 ) 

the exact variation depending on the particular tire 
force-deflection characteristics. Prior to the beginning of 
shock-strut deflection 

Fr=F V[ (z s ) (la) 

since z u =z f . (This relationship is exact when the landing 
gear is vertical and holds very closely when the gear is 
inclined.) 

The shock strut starts to deflect at the time t T when the 
force exerted on the airplane by the shock strut becomes equal 



(b) 

(b) Schematic representation of shock strut, 
airplane) considered in reference 6. 

to the air-pressure preloading force in the strut. At this 
instant the free-body equation for the unsprung mass of the 
landing gear is 

m u z, r +F Vt (z fT )=Pa 0 Aa cos 8+W, (t=U) (2) 

Equation (2) provides the relationship between the ter- 
minal conditions for the first phase of the impact which, in 
conjunction with the solution of the equations of motion for 
the complete system prior to shock-strut deflection, deter- 
mines the time t T when the shock strut begins to deflect and, 
thus, the terminal values of the variables for the first phase of 
the impact. These values then serve as the initial conditions 
for the second phase of the impact. 

After the shock strut begins to deflect, the landing gear 
has two degrees of freedom, since the motions of the landing- 
gear attachment point and the motions of the unsprung 
mass are no longer the same. The equation for the vertical 
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componentof the force transmitted to the airplane by the land- 
ing gear after the shock strut starts to deflect is (see ref. 6) 


_|i| 2 (G4„) 2 **+*« A (a-aJ) _ 005(1 {t - fr) (3) 


where 


cos 8 

. Zf g u 
cos 8 


The equation of motion of the unsprung mass is 

m u 2 u +F V[ (2 u )=F Vi +W u (t^U) (4) 

t 

In equation (3) the first term represents the hydraulic 

force in the shock strut, where the factor ,-rr indicates the 

1*1 

change in sign required between the compression and exten- 
sion strokes. (During the extension stroke of the shock 
strut, because of the action of the rebound check valve or 
“snubber” incorporated in most landing gears, the net orifice 
area will generally be smaller and the orifice discharge 
coefficient C d will be different from the values which apply 
during the compression stroke.) The second term of equa- 
tion (3) expresses the air-compression force in the strut, 
based on a polytropic pressure-volume relationship. In 
equation (4), the force arising from the deflection of the 
tire may be expressed as F Vg (z u )=Tnz u T for the usual types 
of pneumatic tires, where to and r are constants for each 
regime of the tire-deflection process (see ref. 6). 

EQUATIONS OF MOTION OF THE AIRPLANE 

Differential equations of airplane structure. — In the mode- 
superposition approach, the structure is considered to deflect 
in its natural modes of vibration and the total displacement 
of any point in the system is the sum of the displacements 
of the point in all the modes considered. With this approach 
the motions are separated into functions which depend only 
on the space coordinates and functions which depend on the 
time variable. 

In the case of a landing impact the process is discontinuous 
at the instant U when the shock strut begins to deflect. In 
the first phase of the impact the shock strut is effectively 
rigid, so that the motion of the unsprung mass of the land- 
ing gear is essentially the same as the motion of the landing- 
gear attachment point and the force transmitted by the 
landing gear to the airplane is the vectorial sum of the 
ground force due to tire deflection, the inertia reaction of 
the unsprung mass, and the weight of the unsprung mass. 
In the second phase of the impact, the motion of the un- 
sprung mass is not the same as the motion of the landing- 
gear attachment point and the force applied to the airplane 
is governed by the relative motion between the landing-gear 
attachment point and the unsprung mass, as given by 
equation (3). 



(O) 




(a) Coordinates along elastic axis. 

(b) Coordinates at any transverse station. 

(c) Coordinates at landing-gear station. 


Figure 2. — Coordinates for airplane structure. 

The notation employed in the analysis is indicated in 
figure 2. A typical transverse station located at a span- 
wise distance y from the airplane center plane is considered. 
The mass per unit length of span is designated by to. The 
translation of the elastic axis at the station is denoted by w; 
f is the translation of the station mass center; e is the chord- 
wise distance between the station mass center and the elastic 
axis; and <p is the angle of twist of the station. The trans- 
lation of an arbitrary point located at a chordwise distance 
x from the elastic axis is designated by z. The spanwise 
distance from the center plane of the airplane to the landing- 
gear station is indicated by y t . The translation of the 
landing-gear attachment point, or force-application point, is 
designated z f ) the distance between the landing-gear attach- 
ment point and the elastic axis is denoted by x f . 

In the most general case, the expansion of the deflection 
of the structure in terms of its natural coupled modes of 
vibration may be written as 

w(y,t)=^2a*(t)w n (y) (5) 

n-0 

and 

<p(y,t)='^ l a x {t)<p H (y) (6) 

n *=*0 
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where the subscript n denotes the order of any mode, a n is 
the generalized coordinate in the 77 th mode, and w n and <p H 
are the corresponding modal functions for bending and 
torsion, respectively. 

For later use it is convenient to introduce expressions for 
the displacements at other points in the structure. Since 
the translation of the station mass centers is given by 
t=w-\-6<p, it follows that 

■ n — 0 

where the modal function ^ n =w n J r e<p n . The translation of 
any arbitrary point along the chord is given by z=w+x<p; 
therefore, 

z(x,y,t)==i?fi*(t)z u (y) (8) 

n «=»0 


where the modal function z n =w n -\-xtp, ,. The translation of 
the landing-gear attachment point is given by z f =w-\-Xf<p] 
therefore, 

(9) 

n «0 

where the modal amplitude £ n =w n +X/<p n . 

By application of Lagrange’s equation and the orthogo- 
nality relationships between coupled modes, it can be shown 
(see, for example, refs. 8 to 10) that the equation of motion 
for the airplane in the Tith mode may be written as 

M n tf B +M„u* 3 a,= <?» (71=0,1,2, . . .) (10) 


where M n is termed the generalized mass for the 7ith mode 
and Q n is the generalized force, as determined from virtual- 
work considerations. For a continuous system, 


fb/2 fb/2 rb/t 

M a = mw x 2 dy+2 m&p,w n dy+ mK 2 <p x dy 
Jo Jo Jo 

rbn rbn 

= rn$ 2 dy+ IWdy 
Jo Jo 


( 11 ) 


The relationship between Q n and the external forces can 
be determined by application of virtual-work principles. 
By definition, the work done in the nth. mode by the general- 
ized force acting through a virtual displacement of the gen- 
eralized coordinate of the mode is equal to the work done by 
the external forces acting through virtual displacements of 
their points of application in the mode. Thus, the virtual 
work, done by the generalized force in the nth mode is 

SW n =Q n 8a n (12) 


In the case of an airplane during landing, the external 
forces are the distributed lift forces L(y), the distributed 
weights gm(y), and the force F transmitted by the landing 
gear. The virtual work done by these external forces in 
the Tith mode is therefore given by 


8W n =— L8a K z x dy—gJ° m8a x { x dy+F8a x ^ 
~-(f Lz x dy—gj^ mtndy+F £*) (13) 


Equating equations (12) and (13) gives the following 
relationship for Q n : 

<2n=— n Lz n dy~gj o m ?n dy+F^^j 


Therefore, the equation of motion of the structure in the 
77th mode is 

rb/i i'b/2 

M n d n +M n u n 1 a H =—Ft n — Lz x dy+g\ mf, dy 

Jo Jo 

(77=0,1,2,...) (14) 


For the rigid-body mode ( 77 = 0 ), since a> 0 =0 and Zo=fo=£o= L 
equation (14) becomes 


fb/2 

M(fl 0 = — F — I ( L—gm)dy 


In practice the spanwise mass distribution is often approxi- 
mated by breaking up the distribution into- discrete masses 
which are concentrated at a fini te number of stations along 
the span. With this approach equation (11) may be written 
as 


subject to the initial conditions 


o fl (0)=0 

and 

do(0)= W 0 


M x ~^(m i w, i i +2m i e i <p Mi w n( +m i K { 2 <p nt 2 ) 


(11a) 


where the subscript i denotes any spanwise station. 

For 77=0 (rigid-body mode), since Wo=f 0 =l And ?x>=0, 


rbn 

M 0 = I m dy~ 


If the airplane is assumed to be free of oscillations at the 
time of initial contact, 

d B (0)=d B (0)=0 ( 77 ^ 0 ) 

At the instant of initial contact, the airplane may bo 
accelerating, that is, d£ 0 (0) 5 ^ 0 , if the total lift is not exactly 
equal to the total weight. Consideration of the balance of 
forces on the unsprung mass as a free body leads to the 
following equation for the force applied to the airplane by 
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the landing gear at the instant of contact: 

F(0)=W u [^-l] 

Substitution of this relationship into the equation of motion 
for the rigid-body mode gives 


d«(0)= 


rb/i 

ft 7 *— I ( L—gm)dy 


W toi ~L, 


M 0 +m u 


Wt 0 i/g 


where 


so that 


=(l~K L )g 



F(0)=-K L W a 


With this definition of F( 0) and the initial conditions for 
the modes n^O, equation (14) applied to the instant t= 0 
gives 


phases, as previously discussed, the equations of motion for 
each phase must be solved separately, where the initial 
conditions for the second phase are the same as the terminal 
conditions for the first phase. 

Motion prior to beginning of shock-strut deflection. — 
Since the shock strut is effectively rigid in the first phase of 
the impact, the force transmitted by the landing gear to the 
airplane, F in equation (15a), is equal to the ground force 
F v ( z f ) less the inertia reaction of the unsprung mass and 
the weight of the unsprung mass, as may be seen by consider- 
ing the unsprung mass as a free body; thus, 

Ft^t T —F v g (Zf).'\-m u Zf — W u 


so that the motions of the system dining the first phase of 
the impact are governed by the following set of differential 
equations: 


MoOo= — \FyJ>z f )-\-m u z f ^W tot {Kr, — 1)] 


— [F y [ {Zf)-\-Tn'uZf J rW t FK L !)]& 




1 / rbn rbn \ 

a *(°) = M^r 2 v Jo m ^ ndy ~Jo Lzx dyJrKr ^ u ^) ( n ^°) 

This relationship indicates that, in general, a finite static 
deflection in the flexible modes will be present at the time 
of initial contact. At any subsequent time the deflection 
will be equal to this initial static deflection plus an addi- 
tional deflection a„ ( which varies with time; that is, 
<x„=<x n (0) This substitution permits equation (14) to 

be written as 

M H d Ht +M n b3 n 2 a ni = — (F+KlWu)^ (n^O) (15) 

subject to the initial conditions 

a* t (0)=d n( (0)=0 

In the remainder of the report, for the sake of simplicity 
of notation, the subscript t will be dropped, with the under- 
standing that a H represents the time-varying part of the 
displacement of the nth mode, so that equation (15) is 
written as 


M m a m +M m a m 2 a n =—[F Vg (z f ) +ra u ~z f +W u {Kr~lM n J 


where 


Z J S a n?it 

»*=>0 


(16) 


and the mth mode is the highest mode considered. 

The initial conditions for equations (16) are the conditions 
at the instant of initial contact, namely, 

do(0)=0 
«o(0)= Vy o 

and 

a n (0)=a B (0)=0 0) 

As previously indicated, the first phase of the impact 
terminates at the time U when the force in the shock strut 
becomes equal to the air-pressure preload force. The ter- 
minal conditions at this instant, as determined by consider- 
ation of the unsprung mass as a free body, are given by 
equation (2), namely, 


M x d n + M % u n 2 a n = — (F+K L W v )Z a (n^Q) (15a) 


muZf T +Fv t {zf T )='pa 0 A a cos d+W u 


If the external forces are specified solely as functions of 
time, the equations of motion for each mode of the system 
are uncoupled and can be solved individually. However, 
when the external forces depend on the motions of the sys- 
tem, as in the case of the landing-gear force during a landing 
impact, the relationships between the external forces and the 
motions in the modes serve to couple the equations of motion 
so that they must be solved simultaneously. Furthermore, 
in the case of landing impact, since the process has two 


The solution of equations (16) in conjunction with equa- 
tion (2) permits the determination of the time t, when the 
shock strut begins to deflect and the values of the motion 
variables at this instant; these values then serve as the initial 
conditions for the second phase of the impact. 

Motion subsequent to beginning of shock-strut deflec- 
tion. — In the second phase of the impact the force trans- 
mitted by the landing gear, F in equation (15a), is the ver- 
tical component of the shock-strut force F v , as given by 
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equation (3). Thus, the motions of the system during the 
second phase of the impact are governed by the following 
set of differential equations: 

H/o5o= — (FV^-j-W,,) — W lol (K L — 1) 

il/i<ii+AfiWi*aj= (F Vt - 

and 

m u z n J rFy t {z^=F Vt -\-W u 

where 

Fv=F Vt (z f —z U) Zf—Zu) 

as given by equation (3); and 

m 

Zf— y ) ttn£n 
n>=0 

The first m equations of equations (17) represent the mo- 
tions of the airplane structure in its first m modes, whereas 
the last equation of the set is the equation of motion of the 
unsprung mass of the landing gear as previously given by 
equation (4). The initial conditions for equations (17) are 
the terminal conditions for equations (16) as previously 
discussed. In view of the fact that the landing-gear forcing 
term F v , is highly nonlinear, analytical solution of the sys- 
tem of equations (17) does not appear possible, so that it is 
necessary to resort to numerical-integration or analog 
methods. 

SIMPLIFIED SYSTEM CONSIDERED IN NUMERICAL STUDIES 

The preceding section has presented the equations of 
motion for a flexible airplane coupled to a landing gear, 
which permit calculation of the motions of the system 
during a landing impact with consideration of as many modes 
as may be desired. For the study of the effects of inter- 
action between the landing gear and the structure, however, 
it appears that the primary effects of structural flexibility 
on the behavior of the landing gear can be represented by 
considering only the first deflection mode in addition to the 
rigid-body mode. 2 This simplification, which greatly reduces 
the amount of computational work, is felt to be justified for 
the purposes of the present investigation since both theoret- 
ical considerations and experimental data indicate that the 
higher modes should have relatively little effect on the 
landing-gear performance. With this assumption the equa- 

1 In a dynamic analysis, stresses In the structure due to excitation of the higher modes can 
be approximated by calculating the response of such modes, individually, to the forcing 
function determined for the landing gear coupled with the rigid-body and first deflection 
modes. This procedure should be a considerable improvement over the use of the rigid-body 
forcing function as a basis for response calculations in cases where the landing-gear attachment 
points experience appreciable deflections relative to the mass center of the system- 



tions of motion reduce to 

Modo= — [F Vt (z/)+m u z f +W lol (K L — 1)] ") (18a) 

M l d l +M 1 o> 1 i a 1 =— [F Vg (z r )-\-m u z f -\-W u {K L — 1)]£, J (18b) 

and 


Flgda — — (F^+Wu) W tot (B lx, 1) (19ft) 

M 1 d 1 -\-Ad l (jy l 2 a\= — (F V) -f- K L Wv) £ i K (19b) 


where 

and 


m u z u +F Vg (z u )=F Vg +W u 

2/=Oo + Ol5l 

m u z fr +F Vt (z fr )=p ao cos 0+ W u 


0 t=tr ) 


(19c) 


The solution of equations (18) and the determination of 
the conditions at the time t r when the shock strut begins to 
deflect, which serve as the initial conditions for equations 
(19), are treated in appendix A. With these initial condi- 
tions, equations (19) may be solved by numerical integration 
or analog methods. 

From the time-history solutions for the motions of the 
system thus obtained, the accelerations and inertia loads at 
any point in the structure can be calculated from the equa- 
tions presented in appendix B. 


equivalent three-mass system 

It is of interest to note that the equations of motion 
previously given not only represent the distributed 
system of the airplane but can also be used to define equiva- 
lent systems of spring-connected masses, where the number 
of masses above the landing gear is equal to the number of 
modes considered. For the particular case where two modes 
are considered, the equivalent system is one containing three 
masses, one of which is the unsprung mass of the landing 
gear. The use of such a three-mass system provides a rela- 
tively simple means for simulating the primary effects of 
structural flexibility in actual drop tests of landing gears in 
a drop-test jig. 

In the equivalent three-mass system (see fig. 3), m t repre- 
sents the mass to which the landing gear is directly attached 
and m s is the elastically connected mass. The displacement 
of m r relative to its position at the instant of initial contact 
is denoted by z f ; the displacement of m, is designated z„ 
whereas the displacement of the axle or unsprung mass m„ 
is z«. The spring constant of the elastic member is denoted 
by k. Separate lift forces L, and L f will be considered to 
act on the masses m, and m f . 

In order that the three-mass system represents the airplane 
properly, z f , z u , m v> and, of course, the landing-gear charac- 
teristics must be the same for the two systems, so that the 
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Fiqubb 3. — Equivalent three-mass system for representing flexible 
airplane and landing gear. 


landing-gear force is the same; and m r , m„ k, and the applied 
lift forces must be determined from the relationships between 
the equations of motion for the three-mass system and the 
equations of motion for the airplane. 

Consideration of the forces acting on each mass as a free 
body (see fig. 3) leads to the following equations of motion 
for the three-mass system: 

Prior to beginning of shock-strut deflection, 

(m f -\-m u )z f — k(z,— z f )+L f — (W f +W u ) (20a) 

= F v g {zj) 

r (tgtr) 

w 1 2,+(m / +m u )2 / +(£, +.£./) — (W, -f- W/+ TP*) 

Fy (Zf) (20b) 

* 

where 

m u z fT J r F Vg {z fr )^=j) a A a cos 0+TF u (t=U) 

Subsequent to beginning of shock-strut deflection, 
m f Zf— k{z,— z f )-\-Lf— W r =— F Vt ' (21a) 

m/z f -^-in,z s -\-(L,-\-Lf) — ( W, + TP/) = — Fy t *■ (F>tr) (21b) 

™-v.Zu-\-Fy g (zJ)—F Vt -\rWu (21c) 

The problem is to determine the relationships between 

m„ m f , k, L„ and L f for the airplane so that equations (20) 
are equivalent to equations (18) and equations (21) are 
equivalent to equations (19) with the requirement that the 
motions of the landing gear for the three-mass system be 
the same as for the airplane, that is, 

Zf=a o+ffifi 

and that z u be the same in both systems. Since equations 
(19c) and (21c) are identical, they need not be considered 
further in evaluating the unknown constants for the three- 
mass system. 

It is apparent that equations (20a) and (20b) as well as 
equations (21a) and (21b) can be written as 

mjZ f —k{z,—zj)-\rL r —W^—F (22a) 

m f z f +m,z,+ (L,-{-Lf) — (W,+ W/) = — F (22b) 

where 

F=F Vt (z u )-\-m u z u — W u 

and 

2 tI =2 / (tA.tr) 

Zu^Zf (t>t r ) 

Similarly, equations (18a) and (18b) and equations (19a) 
and (19b) can be written as 


■Mo<Eo+ W tot(K L — 1) + W u = — F 


(23a) 
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^a l +KJV u =-F (23b) 

51 51 

Thus, the problem is reduced to determining the constants 
for the three-mass system so as to make equations (22) 
identically equivalent to equations (23) . This may be done 
in any of several different ways. For example, since the 
structure is taken as linear, let 

z,=ax>~\~o^ 

where /3 is a constant to be determined. Substituting for 
z r and z, in equations (22a) and (22b) and eliminating do 
between these equations gives 

ni/(f i— /3) <h+~~ (Si~ 0)( m /+ w*'j)ox+^'/~ 

W,+^(W.-L,)=-F (24) 

whereas subtracting equation (22a) from equation (22b), 
with the same substitutions, gives 

dk>+/kh+~-~ — li)ffi+ _* — <7=0 (25) 

m, m, 

Equation (24) is directly comparable with equation (23b). 
Combining equations (23 a) and (23b) so as to eliminate F 
and to make the coefficient of do equal to unity gives the 
following equation with which equation (25) may be 
directly compared: 

(26) 

In order to evaluate the constants for the three-mass sys- 
tem, each term in equations (24) and (25) is set equal to 
the corresponding term in equations (23b) and (26), re- 
spectively. This procedure gives six simultaneous equa- 
tions, the solution of which yields the following expressions 
for the constants in the three-mass system: 



(27) 

L,+Lf=K L W, ot 

(28) 

m M 0 V 

(29) 

m ' = M l +M 0 £ 1 * 

(30) 

m. Mol! 2 
m f Mi 

(31) 



k-Mm ( Mi+Afo£i *) 

(32) 



(33) 

and 

L f -K L g ( Mi+Mofi 2+ m «) 

(34) 

where 

» Mi 

P Mo|! 

(36) 

and 

mjz f -\- m,z , = MxJo 

(36) 


With the foregoing substitutions, equations (22) are iden- 
tically equivalent to equations (23); thus, the three-mass 
system with the specified values of m„ m f , k, L„ and L f can 
be considered to be equivalent to the airplane in its first two 
modes during both the first and second stages of the im- 
pact. Equations (27) and (28) are required to satisfy the 
equations of motion for the airplane as a rigid body, whereas 
equations (29) to (34) are required for proper representa- 
tion of the airplane in its first flexible mode. With this ap- 
proach the structural properties of the airplane are defined 
by three parameters: the total mass above the landing gear 
Mo, the mass ratio mjm f , and the natural frequency 

The solution of the equations of motion during the first 
phase of the impact and the determination of the conditions 
at the instant of initial shock-strut deflection t T are treated 
in appendix A. With these conditions as initial conditions, 
the equations of motion for the second phase of the impact 
can be solved by numerical-integration or analog methods. 
From the time-history solutions for the motion of the three- 
mass system, the inertia loads and bending moments at any 
point in the airplane structure can be calculated by use of 
the equations in appendix B. 

SOLUTION OP EQUATIONS OF MOTION 

In view of the fact that the equations of motion subsequent 
to time U are highly nonlinear and therefore cannot bo 
solved in closed form, it is necessary to resort to numerical- 
integration or analog methods. Various numerical- 
integration procedures are given in references 11 to 13. 
Appendix A of reference 6 illustrates the application of 
several such methods to the problem of the impact of a 
landing gear attached to a rigid mass. One of these methods, 
which may be termed the “quadratic procedure,” was used 
to obtain those numerical results presented in this report 
which could not be obtained analytically. 

In this procedure, which involves a step-by-step solution 
of the equations of motion, the following difference equations 
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(ref. 11, p. 16) based on a quadratic variation of displace- 
ment over successive finite time intervals are used to replace 
tlie derivatives in the equations of motion: 


and 


Ap 2 At 

^p+i Xp_ i 

(a ty 


where \ p is the value of any variable at the end of the pth 
interval subsequent to the beginning of the process and At 
is the time interval. The difference equations of motion 
obtained by substituting these expressions into the differ- 
ential equations of the system then become essentially extra- 
polation formulas which permit calculation of the displace- 
ments to come from the values of displacement already 
calculated, the whole procedure starting out with the initial 
conditions of the process; that is, the conditions at the 
instant t=t T when the shock strut first begins to deflect. 
With the displacement time histories thus calculated, the 
velocities and accelerations are then determined from the 
foregoing difference equations. 

CALCULATED RESULTS AND DISCUSSION 


CASES CONSIDERED 

In order to investigate the effects of structural flexibility 
on the behavior of the landing gear and the loads in the air- 
frame, several case-history studies are presented which cover 
a rango of airplane mass ratios m,/m f . The calculations are 
based on the structural properties of two large airplanes 
having considerably different mass and flexibility character- 
istics. Airplane A is representative of a four-engine propeller- 
driven World War II bomber having a gross weight of 47,200 
pounds and a natural frequency of vibration in the first 
coupled bending-torsion mode of 3.37 cycles per second. 
The structural characteristics used for airplane B are repre- 
sentative of a present-day swept-wing six-jet-engine bomber 
having a gross weight of 125,000 pounds and a natural fre- 
quency of 1.29 cycles per second in the first coupled bending- 
torsion mode. The landing-gear characteristics used for 
airplane A were based on the manufacturer’s data, whereas 
for airplane B, because information was not available, the 
shock-strut characteristics were chosen so as to yield a 
landing gear which is essentially a scaled-up model of the 
landing gear of airplane A. The pertinent numerical data 
for airplanes A and B are given in tables I and II, respectively; 
the modal functions for the first coupled bending-torsion mode 
are plotted in figure 4. 

Tho main landing gears of airplane A were located in the 
inboard engine nacelles very close to the nodes of the first 
coupled bending-torsion mode; in the case of airplane B the 




(a) Airplane A. 

(b) Airplane B. 

Figure 4. — Modal functions for bending and torsion. 


landing gear is of the bicycle type and is located in the air- 
plane center plane. The position of the landing gear (since 
it determines the value of the modal amplitude fi) in con- 
junction with the values of M 0 and M x governs the value of 
the mass ratio m,/m / for each case. (See eq. (31).) 

In order to represent a broader range of mass and flexibil- 
ity effects, the calculations for each airplane were made for 
four mass ratios corresponding to three arbitrary landing- 
gear positions in addition to the original landing-gear 
location. In practice, of course, a change in landing-gear 
location would probably necessitate a modification of the 
wing structure and result in some change in the modal 
characteristics and, thus, the mass ratio. The main pur- 
pose of the calculations, however, is to indicate the effect of 
mass ratio on the behavior of the system, and the exact 
locations of the landing gear which correspond to the mass 
ratios used are of secondary interest. 

In the calculation of the mass ratio tho landing- 

gear force was assumed to pass through the mass center of 
the landing-gear station. Since the modal characteristics 
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used were for the complete airplane including the unsprung 
mass of the landing gear m it was assumed that the un- 
sprung mass was rigidly connected to the mass m f in the 
equivalent three-mass system, as in the first phase of the 
impact, so that 

wi, 

Mi 

where M 0 , Mi, and £1 include the effects of the unsprung 
mass as part of the airplane mass distribution. The mass 
ratios considered and the corresponding landing-gear locations 
are as follows : 


Airplano A 

Airplane B 

Landing-gear location at— 

Mass ratio, 
vufmf 

Landing-gear location at— 

Mass ratio, 
vu(mf 


a 24 


0.22 


0 


0 

Station 215 

.52 

8tation 420. - 

.85 

8tatlon307___ — 

3.33 

Station 504 __ 

2. SI 




When the landing gear is located at the node of the first 
flexible mode, this mode, of course, is not excited and, since 
higher modes are not considered in the numerical calculations, 
the airplane behaves as though it were a rigid body, its 
motion being governed by equation (23a). As might be 
expected, the farther away the landing gear is from the nodes, 
the larger is the effective flexibility of the system and, thus, 
the mass ratio. 

In the calculation of the time histories of the motions of 
the system, the lift force was assumed to be constant dining 
the impact and equal to the total weight of the airplano, 
that is, K l = 1. This assumption corresponds to the con- 
dition that 

L,=W. 

and 

L f =W f +W u 

in the equivalent three-mass system. 


TABLE I.— CHARACTERISTICS OF AIRPLANE A 


TABLE II— CHARACTERISTICS OF AIRPLANE B 


(a) Structure 


(a) Structure 


[Data taken from ref. 8] 


[Unpublished data] 


Station, 

in. 

m { , 

lb-sec 2 

in. 

I tail 

lb-in.-sec* 

oc, in. 

w,< 

vu 

0 

28. 5 


0 

-0. 078 

0 

133 

16. 3 

85, 234 

-39. 26 

-. 031 

-. 00084 

217 

6. 27 

1, 288 

0 

-. 047 

— . 0016 

307 

9. 15 

61, 717 

-62. 19 

. 164 

00183 

428 

. 974 

536 

0 

. 374 

00185 

548 

. 686 

287 

0 

. 670 

-. 00187 

638 

. 163 

34. 1 

0 

. 936 

-. 00188 


, , lb-sec 1 

iuo, — : , 

’ in. 

„ r lb-sec 2 

' in 


i. cps. 


61. 033 

L 607 
3. 365 


(b) Shock strut 


Station, 

in. 

m, 

lb-sec 2 

in. 

I*ai) 

lb-in.-sec 2 

«/, in. 



0 

109. 534 

4, 475, 280 

200. 37 

-0. 0686 

-0. 000170 

84 

A 695 

3, 046 

-A 65 

-. 0579 

-. 000187 

MB 1 1 

A 920 

19, 490 

-24. 20 

-. 0350 

-. 000204 


22. 177 

278, 942 

-101. 22 

. 0037 

000231 


2. 560 

2, 161 

2. 44 

. 090 

-. 000272 

420 

2. 557 

1, 988 

2. 00 

. 1842 

-. 000322 

504 

L 773 

1, 136 

. 92 

. 3263 

-. 000379 

588 

3. 269 

2,474 

-14. 79 

. 4772 

-. 000436 

672 

8. 628 

8,439 

-26. 88 

. 6309 

-. 000482 

756 

1. 144 

500 

. 60 

. 8181 

-. 000614 

840 

. 520 

186 

6. 48 

1. 000 

000626 


Mo, 


lb- sec 8 
in. 


161. 776 


Mi, 


lb-sec 2 


6. 9006 


[Manufacturer’s data] 

A*, sq ft 0. 163 

A a , sq ft 0. 214 

A„ sqft 0. 00173 

v 0 , cu ft 0. 2597 

p a<j , lb/sq ft 30, 528 

P, slugs/cu ft L 626 

(c) Unsprung mass 

[Manufacturer’s data] 

IF., lb 700 

Tires (one per landing gear) 56-inch smooth contour 

Tire pressure, lb/sq in. 70 

m, lb/ft 86, 309 

r - 1 22 


fi, cps 1. 20 

(b) Shook strut 

[Values estimated from generalized curves of ref. 6] 

2 T&dAnp ’ slug8 / ft 17, 900 

A a , sq ft 0. 686 

v 0 , cu ft 0. 7096 

Pa a , lb/sq ft 30, 628 

(c) Unsprung Mass 
[Manufacturer’s data] 

TV., lb 2,300 

Tires (two per landing gear) 66 X 16 

m per landing gear, lb/ft 280, 180 

r 1.21 
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On the basis of the calculations in reference 6, the shock- 
strut orifice discharge coefficient Ca was assumed as 0.9 
and the polytropic exponent n for the air-compression 
process was taken as 1.12. 

EFFECT OF INTER ACTION ON BEHAVIOR OF SYSTEM 

Time-history solutions for the motions of the system 
during impact at an initial vertical velocity of 10 feet per 
second have been made for the eight configurations previously 
mentioned. Figures 5 to 8 show the variation during im- 
pact of the more important quantities, such as the landing- 
gear force F, the responses <h/g, ajg, z f /g, zjg, the landing- 
gear-motion variables, and the accelerations at the mass 
centers of several stations along the span. Comparison of 
the calculated results for the flexible cases with those for 
the airplane as a rigid body (or landing gear at nodes, 
7?!,/%= 0) indicates that the interaction between the flexible 
structure and the landing gear can result in an appreciable 
reduction in the applied landing-gear force (and thus, the 
nodal acceleration), the largest reductions occurring at the 
highest mass ratios. Furthermore, the reductions in landing- 
gear force at the higher mass ratios are greater for airplane 
B, becauso of its lower natural frequency, than for airplane A. 

Consideration of the calculated time histories of the mo- 
tion of the landing gear indicates how the interaction be- 
tween the flexible structure and the landing gear affects the 
loads produced in the landing gear. Because of the flexi- 
bility of the structure, the landing-gear attachment point 
deflects upward relative to the nodes, or instantaneous 
center of mass of the system, as the applied force builds 
up and the deceleration of the landing-gear attachment 
point is greater than in the case of the rigid airplane. Thus, 
the downward velocity of the shock-strut outer cylinder is 
more rapidly dissipated and the displacement of the outer 
cylinder is smaller throughout most of the impact. The 
tire deflection is also smaller; however, because of the high 
stiffness of the tire, the decrease in tire deflection is smaller 
than the decrease in outer-cylinder displacement. The 
net result is a reduction in strut stroke during that part of 
the impact when the maximum force occurs and an accom- 
panying reduction in the strut telescoping velocity. Since 
the maximum landing-gear force is primarily due to the 
hydraulic resistance in the strut (because the strut stroke, 
and thus the air-compression force, is generally small at 
the time of maximum telescoping velocity), the decrease in 
telescoping velocity results in a decrease in shock-strut 
force. 

In the case of airplane A with landing gear at station 
307, the effect of interaction is a marked change in the shape, 
as well as in the magnitude, of the time histories. Because 
of the superimposed vibrations of the structure, the shock- 
strut telescoping velocity (see fig. 5) has acquired an os- 


cillatory character with two peaks of the same amplitude. 
However, since the second telescoping-velocity peak occurs 
when the stroke is large, the superposition of the high air- 
compression force on the hydraulic-force results in a total- 
force time history the second peak of which is much higher 
than the first (see Force-time curves, fig. 5) and which is 
also higher than might be expected from the results for the 
smaller mass ratios, which have a considerably different 
appearance. In the case of airplane B, because of the 
lower natural frequency, this double-peaked characteristic 
does not appear even for the largest mass ratio, all mass 
ratios yielding time histories similar in shape, the maximum 
force decreasing in a regular manner with increasing mass 
ratio. 

The extent to which the first flexible modes of airplanes 
A and B are excited by the impacts may be observed by 
examining the time histories of a h d u and a x . As may be 
expected, the higher the mass ratio, the greater is the degree 
of excitation. 

From the calculated values of a 0 /g and djg or zjg and zjg, 
the acceleration at any point along the span may be com- 
puted by means of the equations in appendix B. Figures 6 
and 8 show time histories of the acceleration at the mass 
centers of several stations for each of the landing-gear loca- 
tions considered. Because of the combined effects of the 
changes in the landing-gear forcing function and in the degree 
of excitation of the flexible modes, a given change in landing- 
gear location may result in an increase in acceleration at 
some stations and a reduction in acceleration at other sta- 
tions. 

Figures 5 and 7 also show time histories of the acceleration 
zjg which would be experienced by the elastically connected 
mass m, in the equivalent three-mass system, as in a drop 
test. The reduction in acceleration with increasing mass 
ratio is evident. As previously indicated, if such a drop test 
were made, the measured accelerations zjg and zjg could be 
used to calculate the accelerations and stresses that would 
result at any point in the corresponding airplane structure 
by means of the equations presented in appendix B. 

Figure 9 (a) presents a summary graph showing the effects 
of structural flexibility and interaction on the maximum 
landing-gear force for the various configurations considered. 
As previously indicated, the reductions in landing-gear force 
are greater for airplane B than for airplane A because of the 
lower natural frequency of airplane B. For the range of mass 
ratios representative of existing and proposed large air- 
planes, for example, values up to about 0.5, reductions in 
landing-gear force up to between 15 and 20 percent may be 
possible. Along the same lines, figure 9 (b) shows the effects 
of interaction on the acceleration response of the landing-gear 
attachment point and on the acceleration of the elastically 
connected mass in the equivalent three-mass system. 


430875—57 41 



Landing-geor force, F , lb 



0 J 230 J 230 J 2 3 
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Figure 5. — Time histories of motions of airplane A during impact. 
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(a) Landing-gear force. 


Moss ratio, m s /rrif 


(b) Accelerations of landing-gear attachment point and elastically connected mass. 
Figure 9. — Effects of interaction. 
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EFFECTS OF NEGLECTING INTERACTION IN THE 
CALCULATION OF DYNAMIC LOADS 

In the usual procedures of dynamic analysis of landing 
loads it is customary to neglect the effects of interaction on 
the landing-gear forcing function and to determine the dy- 
namic loads in the structure by calculating the response of 
the structure to the forcing function which would be obtained 
if the airplane were a rigid body, this rigid-body forcing 
function being either calculated or, more frequently, de- 
termined on the basis of drop tests of the landing gear with a 
rigid mass. In practice, either the actual rigid-body forc- 
ing function or some simplified analytical approximation of 
it (see, for example, fig. 10) is used. 

In order to evaluate the errors introduced by neglect of 
interaction effects, the root bending moments and shears 
determined from the interaction solutions for airplanes A 
and B are compared in figures 11 and 12 with those deter- 
mined by calculating the response of the various configura- 
tions to the rigid-body forcing functions previously presented 
and to simple analytical approximations to the rigid-body 
forcing functions. These bending moments and shears are 
* total values due to both inertia and aerodynamic forces, the 
latter being included to permit comparison with the steady- 
llight values. For reference purposes, figures 11 and 12 
also show the root bending moments and shears which would 
be experienced by a completely rigid airplane. 


The calculation of the response of systems with two de- 
grees of freedom to predescribed forcing functions is treated 
in appendix C. The response of the various configurations 
to the rigid-body forcing function was calculated by applica- 
tion of the numerical-integration procedure previously 
described, whereas the response to the analytical forcing 
functions was obtained in closed form. The rigid-body 
forcing functions for airplanes A and B and their approxima- 
tions are shown in figure 10. In the case of airplane A, 
the rigid-body forcing function was approximated by a pulse 
composed of sine and cosine segments; for airplane B, a 
simple sine pulse was used. The equations for calculating 
the inertia moments and shears from the response of the 
system are given in appendix C; simplified expressions for 
calculating the moments and shears due to the aerodynamic 
forces are given in appendix D. 

From figures 11 and 12 it can be seen that the bending 
moments and shears calculated from the response to the 
rigid-body forcing function are larger than those determined 
from the interaction solutions, the differences being greater 
for the higher mass ratios where the effects of interaction 
result in a greater reduction in the magnitude of the landing- 
gear forcing function. From these particular examples, it 
appears that neglect of the effects of interaction on the 
landing-gear forcing function can lead to overconservatism in 
design not only of the landing gear but also of the structure, 



Time after contact, sec Time after contact, sec 

(a) Airplane A. 0=12.08 radians/second; Oi=8.27 radians/second. (b) Airplane B. Q=12.57 radians/second. 

Figure 10. — Rigid-body forcing functions and simple analytical approximations. 
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(a) Bending moments. 

Figtjeb 11. — Dynamic loads in airplane A. 

particularly for very flexible configurations with high mass 
ratios. As might he expected, there was relatively little 
difference in the loads calculated from the response to the 
analytical approximations and from the response to the 
rigid-body forcing function. 

It is of interest to note that in the case of airplane A the 
loads calculated from the interaction solutions are greater 
than those calculated for the completely rigid airplane, 
whereas, for airplane B, the converse is true. This result for 
airplane B is due to two factors : (a) the dynamic amplifica- 
tion factor is less than unity because of the relatively large 
natural period of the airplane compared with the duration of 
the impact pulse (£j/£*«D.3), and (b) there is considerable 
reduction in the magnitude of the landing-gear force because 
of the effects of interaction. In the case of airplane A, the 
natural period is of about the same duration as the impact 
pulse (£«/£»*» 1.1) so that the dynamic magnification factor is 
considerably greater than unity and more than overcomes the 
effect of the reduction in landing-gear force. 

From the preceding results, it can be seen that the effects 
of structural flexibility are twofold; namely, (a) a change in 



Figure 11. — Concluded. 

the magnitude of the applied landing-gear force duo to inter- 
action, the amount depending on the natural frequency of 
the structure, the mass ratio m,/m/, and the landing-gear 
characteristics, and (b) either dynamic amplification or 
attenuation of the loads in the structure compared with 
those for a rigid body, depending largely on the ratio of the 
duration of the impact pulse to the natural period of the 
structure. In the particular examples considered, the 
landing-gear force was reduced by the effects of interaction ; 
it is conceivable, however, that, for some combinations of 
landing-gear and airplane characteristics, perhaps when the 
natural period of the structure is smaller than the duration 
of the impact pulse and the mass ratio is large, interaction 
might result in an increase in the maximum landing-gear 
force over that for a rigid airplane because of the superposi- 
tion of oscillations of the landing-gear attachment point on 
the motions of the shock strut. Such an unfavorable effect 
of structural flexibility of the applied force was indicated for 
certain cases of seaplane impact in reference 7. 

In view of the foregoing observations it would appear 
worthwhile to consider the effects of interaction in dynamic 
analyses of landing loads when the landing gear is located 
at points in the airplane that experience appreciable defloe- 
tions relative to the mass center of the system. 
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(a) Bending moments. 

Figube 12. — Dynamio loads in airplane B. 




Time after contact, sec 


(b) Shears. 

Figube 12. — Concluded. 


CONCLUSIONS 


The effects of interaction between a landing gear and a 
flexible airplane structure on the behavior of the landing 
gear and the loads in the structure have been studied bj 
treating the equations of motion of the airplane and the 
landing .gear as a coupled system. The landing gear is 
considered to have nonlinear characteristics typical of 
conventional gears, namely, velocity-squared damping, 
polytropic air-compression springing, and exponential tire 
force-deflection characteristics. For the case where only 
two modes of the structure are considered, an equivalent 
tliree-mass system is derived for representing the airplane 
and landing-gear combination, which may be used to simu- 
late the effects of structural flexibility in jig drop tests of 
landing gears. 

As examples to illustrate the effects of interaction, numeri- 
cal calculations, based on the structural properties of two 
largo airplanes having considerably different mass and flexi- 
bility characteristics, are presented. In order to cover a 
range of parameters, the landing gear of each airplane was 
assumed to be located at three arbitrary spanwise positions 
in addition to its original location. For the particular cases 
considered, it was found that 

1. The effects of interaction can result in appreciable 
reductions in the magnitude of the landing-gear force, 
particularly when the flexibility of the airplane structure is 

435875—07 42 


large and the natural frequency of the structure is small. 

2. Neglect of interaction effects, that is, the use of the 
landing-gear forcing function for a rigid airplane in a dynamic 
analysis of a flexible airplane, can lead to the calculation of 
excessive loads in the airplane structure. 

3. In the case of one of the airplanes, the structural loads 
calculated from the interaction solutions are greater than 
those for a completely rigid airplane treatment (rigid struc- 
ture subjected to rigid-body forcing function) because of the 
fact that the effects of dynamic magnification more than 
overcome the reduction in landing-gear force due to inter- 
action. In the case of the second airplane, because of the 
relatively large natural period of the structure in comparison 
with the duration of the impact pulse, the dynamic magnifica- 
tion factor is appreciably less than unity. This effect, 
coupled with the reductions in landing-gear force due to 
interaction, results in structural loads that are less than 
those for a rigid airplane. It thus appears desirable to con- 
sider the effects of interaction in dynamic analyses of landing 
loads for large airplanes, particularly when the landing-gear 
attachment points experience large deflections relative to the 
mass center of the airplane. 

Langley Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Langley Field, Va., May 5, 1955. 



APPENDIX A 

CONDITIONS AT BEGINNING OF SHOCK-STRUT MOTION 

Since the shock strut does not begin to deflect until the preloading force imposed by the internal air pressure is overcome by 
the inertia forces, the shock strut is essentially rigid during the interval between the instant of initial contact with the ground 
and the beginning of shock-strut motion at some time t=tr . During this interval, since the deflection of the tire is essentially 
the same as the displacement of the landing-gear attachment point, the system used in the numerical calculations to represent 
the airplane and landing-gear combination has only two degrees of freedom, namely, the rigid-body or zero-mode displacement 
and the deflection in the first flexible mode, the higher modes being neglected. The purpose of this appendix is to consider 
the motions of the system prior to the beginning of shock-strut deflection in order to determine the conditions which exist 
at the instant the shock strut first begins to move; these motions then serve as the initial conditions for the equations of motion 
of the system during the main part of the impact. For this purpose it may be reasonably assumed that the tire force- 
deflection relationship is linear for the relatively small range of deflection prior to the beginning of shock-strut motion and 
that, therefore, F V[ {z f ) —m'z/. In order to avoid a step jump in the time-history solution at the time U, the constant m' 
should be determined so that 

m’ Zf T =mz f * (Al) 

DISTRIBUTED SYSTEM 

Prior to time U the equations of motion for the airplane and landing gear are given by equations (18) with initial conditions: 

%(0) =Oo(0) =Oi (0) = 0 

2 /(0)=do(0)=Vr o 

di(0)=0 


Since (h= Z - y 0<i > equations (18) can be written as 

<i 


Af, 

li 


y 2 


M {1 d 0 =—m'z f —m, u z f —W u>t {K l ,— l) 

~y~ ao= (pF+ m «) Zf+WJJC L — 1) 


The exact solution of equations (A2) can be shown to be 


where 


Vv ° sia At ~^~TT' sin Bt )+ D ^aT cos At ~^B^ cos Bt+O ( 35 -jg*)]}- 


A= 

B= 


^ E-^W-4JF 




„ MjtVCMo+mJ 
0 6 

D=(K L —l)g 


F-- 


m'Miuii 

Q 


(A2a) 

(A2b) 


(A3) 


040 


(7=iHi( A/o-{- wi„) h^oA^oli 1 
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By successive differentiation of equation (A3), the higher derivatives of z f (t) are found to be 

Vy^A'-O cos At—{B 3 —G) cos Bt]+D [(^=^) sin sin (A4) 

~ z ' {t)= A?=w{ V Vo [B(B*-C) sin Bt-A(A*-C) sin Ai)+D[(B*-C r ) cos Bt-{A*-C) cos At] | (A5) 

Mt )= cos £f-A 2 (A*-<?) cos At]+Z>[A(A 2 -<?) sin At—B{B-—C) sin Bt] (A6) 

At the time t T , the equation of motion of the unsprung mass of the landing gear as a free body is given by equation (2) 
which, with F Vs =m'z f , may be written as 

m u z, T +m'z fr =p a A a cos 6+W u (A7) 

Substituting for z /r and z /r in equation (A7) gives a relationship between t T and m' : 

^ 2 — ^ 2 -^ Vv 0 sin At— — ~ - (m'~ m.B 2 ) sin BtJ+Z) (m'—in^A 3 ) cos At— 

-gr g (m'~m u B 2 ) cos Bt+m'C gi)] =Pa e A. cos 0+W u (A8) 

Because equation (A8) is transcendental in both t T and m’ (to' being involved in the constants A and B), in order to obtain 
an explicit solution for U or m', some approximation to the trigonometric terms is necessary, the ordor of the approximation 
depending on the accuracy required. For the determination of U and m' it will generally be sufficient to assume first-order 
approximations for the trigonometric terms where only the first terms of then- series expansions are used. With these approxi- 
mations the solution of equation (A8) for U is 

OfpcAa COS 8+K l Wu) 

tr ~ m'M l MoV Vo (A9) 

As indicated previously, to' cannot be chosen arbitrarily but must be determined in accordance with equation (Al), which 
may be written as 

m' =mz fT r ~ 1 

The first-order approximation for z fr , obtained from equation (A3), is 

Z fr= V V 0 tr (AtO) 

With these substitutions equation (A9) may be written as 


and the equation for m' becomes 



(All) 


(M2) 
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The first-order approximations for the derivatives of z s at time U, from equations (A4) to (A6), are 


Zf T —Vy o —Dtr 

(A13) 

Z,= Vr o [C- (A*+B>)]U-D 

(A14) 

and 


z’/ r = (Vy-m [C- (A 2 +5 2 )] 

(A15) 

With the values of U and in' calculated from equations (Al 1) and (A12) , the v alues of Z;= z u ^, 

z s r ~ i u T , and Zf r —z Ur can bo 


calculated from equations (AlO), (A13), and (A14), respectively. These values provide two- thirds of the initial conditions 
for the process subsequent to the beginning of shock-strut deflection (eq. 19). The remaining initial conditions, for example, 
Oo r , do r , and &o t , can be obtained by manipulation of the differential equations (A2) . From equation (A2b) it can bo seen that 


By differentiation, 


a<v= 


MW 


Oo,- 


[(f? +m “) if * + (?rr +m ') Zt ~i? <+ w «( k l-v] 

Kf +"■) -f <] 


MW 


(A16) 

(A17) 


where, from equation (A2a), 


Differentiating equation (A2a) gives 




m'z fr + m u z, T +W l0 ,(K L — 1) 

Wo 

VatAa cos d+W lot {K h ~\)- i rW u 

Wo , 


«0 = 


m'z fT +m u z, T 

Mo 


(A 18) 


(A19) 


The substitution of equations (A18) and (A19) and the initial conditions previously determined (z /r , z/ rt Z/ r , and z/ T ) 
into equations (A16) to (A18) provides the remaining initial conditions for the second phase of the impact. 

EQUIVALENT THREE- MASS SYSTEM 

The equations of jnotion for the equivalent three-mass system prior to the time t r are equations (20) with initial conditions 

2 /( 0 ) = z ,( 0)=0 

and 

z f {0) = z,(0') = V Vo 


Since it has been shown that equations (20) are identically equivalent to equations (18) for the distributed system when 
the relationships between the constants of the two systems are as defined by equations (27) to (34), it follows that equations 
(A3) to (A15) are equally valid for the three-mass system when the constants are redefined in accordance with equations (27) 
to (34). The redefined constants, in terms of the properties of the three-mass system, may be written as 

c= mW(M o+mj 

Moiuif+mu) 


D=(K L -l)g 
m' , 


E= 


m f +m K 


C 


F= 


m f m'coi 

'Mo(m f -\-m v ) 
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wliero 

and 


M 0 =m / +m l , 



fcMp 

m f m. 


The equations for U and m' , equations (All) and (A12), become 


and 


where 


f L cos e+K L W u )j' 

(A20) 

m'=m llT j^l?(j3 ao A a cos 

(A21) 

jj. m,+m u 



The values of U and m' given by these equations permit the calculation of Z/ r = 2 « r , z f = z « r , and z f = z „ r by means of equa- 
tions (AlO), (A13), and (A14). The remaining initial conditions for the second phase of the impact, z , and its derivatives at 
the time t T , can be obtained by manipulation of the differential equations (20). Solving equation (20a) for z, at time U gives 


Z -=k 


(m/+m Jz /r +(i-|-wi , )z/ r +A/+it/ 
Differentiating equation (20a) and substituting F Vf (z / )=m'z / gives 

i »r=\ ^(m / +mu)z /r +(fc-m , )z /r J 

An expression for z, T can easily be obtained from equation (20b) as follows: 


-w' 


(A22) 


(A23) 


Zj T =' ~~ (m / +m K )z /T +m'z /r +(L,+Z / )—(W,+W / +W u ) 

TTlrf __ 


(A24) 


Equations (A22) to (A24), in conjunction with the values of z / t , z /t , and Zf r previously determined, supply all the initial 
conditions for the second phase of the impact of the equivalent three-mass system. 



APPENDIX B 


DYNAMIC LOADS IN AIRPLANE STRUCTURE 


The equations of motion of the airplane have been pre- 
viously presented in several forms so that solutions for the 
motions of the structure can be obtained in terms of the 
variables a fl and a u and z f , or z f and z,. The purpose of 
this appendix is to present equations from which the accel- 
erations, bending moments, and shears at any point on the 
airplane structure can be calculated once the time-history 
solutions for the basic variables have been obtained. 

< 

ACCELERATION 

At any point. — The absolute displacement at any point on 
the structure (see fig. 2) is 


Since 

and 


Z=W+X<p 


w=a 0 +a 1 w l 




where w x and <p t are the modal functions for bending and 
torsion, respectively, 


and 

Since 


z—a 0 -\-ai{wi-\-x<px) 


z=df 0 4-d 1 (wi+j^ 1 ) 
Zf &0 


(Bl) 


the acceleration at any point may also be written as 


- r I / “ + M 

z=do~r\^' — do) ~ ^ 

Since, from equation (36), 


(B2) 


Along station mass centers. — At the mass center of any 
station the displacement is designated f and a>=« so that 
equation (Bl) becomes 

(B7) 

where f, is the modal function for the station mass center 
and is equal to Wi+e^. Equation (B2) becomes 

f—fio +(2/ — <Jo)|^ (B8) 

Cl 

Equation (B3) becomes 

l=^^mfZf+ Tn »z,+m,{z r —z,)j^ (B9) 

BENDING MOMENTS 

Outboard of landing gear. — The bending moment at any 
spanwise station y } outboard of the landing-gear station 
y t is readily determined by summing up the inertia moments 
produced by the accelerations of the mass centers of all 
stations i between station y } and the tip. Thus, 

m fUVi-yj) (BIO) 

Inboard of landing gear. — The bending moment at any 
spanwise station y } inboard of the landing-gear station y, is 
equal to the sum of the inertia moments produced by the 
accelerations of the mass centers of all stations i between 
station y } and the tip plus the moment produced by tho 
landing-gear force. Thus, 

™ f Uvt-yi)+F(Vt—vi) (Bl 1) 

where 

F= - [Modo+ W t0 ,(K L —l) + W u ] 


, mfZ f +m,z, 

do= M 0 

the acceleration can also be written as 

* = W 0 [] m / S /+ ^z.+m^Zf—z,) (B3) 

Along elastic axis. — At the elastic axis, the displacement 
is designated w and 2=0 so that equation (Bl) becomes 
simply 

(B4) 

Equation (B2) becomes 

w=d(,-\-{Zf — 6o)~z~ (B5) 

Cl 

Equation (B3) becomes 

w= m fZf+ + m * (%— z.) fl (B6) 


SHEARS 

Outboard of landing gear. — The vertical shear at any span- 
wise station y } outboard of the landing-gear station y g is 
sim ply the s um of the inertia reactions due to the accelera- 
tions of the mass centers of all stations i between station 
y t and the tip. Thus, 

(B12) 

Inboard of landing gear. — The vertical shear at any span- 
wise station y } inboard of the landing-gear station y s is tho 
sum of the inertia reactions due to the accelerations of tho 
mass centers of all stations i between station y; and the tip, 
plus the landing-gear force. Thus, 

= m^-lMcSo+W^ (K l - 1)+W u ] (B13) 

i-Jl 
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APPENDIX C 

RESPONSE TO GIVEN FORCING FUNCTIONS 

la this appendix equations are presented for the acceleration response of the airplane structure to predetermined forcing 
functions applied by the landing gear. The cases considered are the arbitrary forcing function, the sine pulse, and a pulse 
made up of sine and cosine segments. For the particular case where the landing-gear forcing function can be represented bv a 
single sine pulse, 

F(t)=F max sin Qt 

where SI is the circular frequency of the applied sine pulse and is expressed by 



where T is the time to reach F max . 

If the forcing pulse is not symmetrical in time about its maximum value, it may be represented by a combined pulse con- 
sisting of a sine function up to the time T and a cosine function subsequent to the time T. This latter function may be written 
as 

F{t , )=F nax cos 0 x i' (t'^0) 

where 

t'=t—T 


and S2i is the circular frequency of the cosine pulse; the initial conditions are the same as the conditions at the time t=T deter- 
mined from the response to the sine-function segment of the pulse. 

The solutions are presented for the distributed system of the airplane (sketch a) and for the equivalent concentrated-mass 
system (sketch b) 



DISTRIBUTED SYSTEM 


The acceleration response of the rigid body or zero mode is immediately evident from the equation of motion^for n— 0, 
namely, 

F(t)-\-W, 0 ,(K L — l)-j-TT M 
Mo 


&o=- 


The response of the deflection modes follows. 

Arbitrary forcing function. — When the landing-gear forcing function is predetermined and arbitrary, the equation of 
motion for the nth mode (eq. (15a)) can be written as 

M n d K +M n ^a n = - [F{t) +K l W u ] (7i7*0) (Cl) 


where F(t) is an arbitrary function of time and a n is the generalized coordinate of the 77th mode. 

The general solution of equation (Cl) may be written as 

ffi*(t)=-yl=- r>(r)sin 0 > a {t-r)dr+ K J^\ n (cos co,t-l)+a„(0) cos co„H-— sin o> n t 
■02*“* Jo n2*“» “* 

The acceleration response is obtained by double differentiating equation (C2) as follows: 

[F ( t)~ Un j F ( r) sin (^{t—^dr+KiWv cos co,tJ— a„(0)w B s cos w n t— d„(0)co* sin co.t 


(C2) 


(C3) 
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Equations (C2) and (C3) are general solutions to equation (Cl) and thus represent the response of any mode to an arbitrary 
forcing function F(t ) . In the present study of landing impact, the initial conditions are 


and 


®»( 0)=0 


d *( 0)=0 


Sine-pulse forcing function. — For the particular case where the forcing function is a sine pulse, the acceleration response, 
as determined from equation (C3), is 

d.(t)=F ~ 2 (ft sin sin Q,t )— sin £ J — +a n (0) o>, s J cos co M t—a n ( 0) co n sin o > K t (C4) 

where, again, a„(0)=0 and d n (0)=0. 

Half-sine — half-cosine pulse. — In this case the response up to time T is given by equation (C4). Subsequent to time 
T the acceleration response, determined from equation (C3), may be written as 

&n(t')=F max ^~- Wni '~ cos ftiO~cos (0) w.’J cos u n t'—a K (0) sin u n t' (05) 

where 

t'=(t—T)^ 1 


a*(0)=a nT 
d»(0 )=a„ T 

EQUIVALENT CONCENTRATED-MASS SYSTEM 

The equations of motion for the concentrated-mass system subject to an arbitrary forcing function are (see eqs. (22)) 


m f z f — k(z,— z/)+L f — W f =—F(t) 
m / 2 / +m.3,+ (L,+L.) - (W,+ W f ) = -F(t) 


Introducing the new variable 


u=z,—z r 

permits the combination of equations (C6) into a single equation in one variable : 

m f u-\-k u==F(t)+J 


where 


J=L,-W,-%(L,-W,) 


The solution of equation (C7), by analogy with equation (Cl), can be written as 

u(£)=-^- J Fir) sin c^(i— (1— cos oji£)-|-u( 0) cos sin co t t 




0>1 


(Cfl' 


(07) 


( 08 ) 


Mo 


where 
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By substituting u(t) for z.—z s in equations (C6) and combining, the following equations for the responses z, and z t can 
be obtained: 

z,(t)= — — f — f F(r) sin wi(t—T)dT-j — — 2 (1— cos coxO+^O) cos sin +— — — (C9) 

v ’ m, {jrifWi Jo m, w? ’ ' ^ J to, K 1 

and 

zJt)=— -f — [F(t) + (L r —W f )\ -\-k \— — f F(t) sin i^{t~ r)dr-\ — — 2 (1— cos co]t)+u(0) cos cojt-f sin coitl"V (CIO) 
TO/ h YJUjWi J 0 VljWi Wl _| J 


In equations (C9) and (CIO), «(0)=«(0) = 0 for the present application to landing impact. 

Sine-pulse forcing function. — For the case where the forcing term is a sine pulse, equations (C9) and (CIO) become 


2.(0 = 


k_ Fmar(Q 8U1 Wit — Wi sin tit) j J. ( 

~~ m, L 


and 


m=— 

m r X. L 71 


rrijw x(fl 2 — aii *) 

sin oit— coi sin 0 1) 
wifOyiiQ ? — c^ 2 ) 


rrijwi 


(1 — cos wjO+^O) cos coxt-j-^^ sin o^t - 

“i 


W-L. 

TO, 


-sin 12 1 


] - {Lf -w f y 


Jk{ 1 — cos tod) 
m /<!>!* 


[m 


,u( 0) - 
cos coit-j — — sin coxt 

COx 


]} 


(Cll) 


(Cl 2) 


where, again, «(0)=ti(0)=0. 

Half-sine — half-cosine pulse. — The response up to time T is given by equations (Cll) and (Cl2). Subsequent to time 
T, the responses (eqs. (C9) and (ClO)) become 


and 


k f Fnar (COS Wit '— COS at') . J 
m, _ to/12x j — &> i 2 ) mw* 


(1— cos wxfO+^O) cos coxt'-f-^^ sin coxt' - 

“i 


W-L, 

TO. 


(C13) 


Zf(t 


'HK'-O 


k (cos wit ' — cos ait') 


to/Qx 2 — o>x 2 ) 


-cos Qit 




Jk{\ — cos wit')., r 


m,jwi 


j^u(0) cos coxt+^^ sin wit 


Wi 


']} 


(C14) 


where 


t'=t-T=0 


u(Q )=Ur 
u{Q)=u t 



APPENDIX D 

AERODYNAMIC AND WEIGHT MOMENTS AND SHEARS 


In appendix B equations were presented for the bending moments and shears due to the combination of the inertia forces 
arising from the accelerations of the masses distributed along the span and the landing-gear force. In the calculation of 
the total moments and shears, however, consideration must be given to the aerodynamic lift and weight forces. This 
appendix presents equations for estimating these aerodynamic and weight moments and shears which, although only first 
approximations, are considered sufficiently accurate for the purposes of the present study. 

If it is assumed that the lift coefficient is constant along the span and equal to the average lift coefficient of the wing 

Cl, the lift force at any station y t is equal to C L ^ V z Mi where A t is the area assumed to he concentrated at the station. 

The moment at any station y } due to the lift and weight forces at each station i outboard of station y } is 

n tlD ‘ up 

. M ttm =C L ^V L i 2 3 4 5 6 'Z^A t (y t —y J )—gZ!,m i (y i ~y J ) (Dl) 

1 f=j {•=) 

If unsteady-state lift effects are neglected, the instantaneous lift coefficient is related to the lift coefficient at the instant of 
initial contact by the expression 


Cl — Cl 0 ~\~ C La (y — %>) 


=CL 0 +C La 




Inasmuch as the total lift at the instant of contact is KlW,,,,, 

Cr- 


K L W tot 


so that 


a PW A 

Similarly, the shear at any station y } is 


(«o -Vy 0 ) | Pi] g My t - yi )~g £ mtfa-y,) 


J^+ C LX a °- V Vo) | Vl~\ fj^i-g § m t 


(D2) 


(D3) 
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